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We consider a long-range Ising antiferromagnet put in a transverse field (LRTIAF) with disorder. 
We have obtained the phase diagrams for both the classical and quantum case. For the pure case 
O ' applying quantum Monte Carlo method, we study the variation of order parameter (spin correlation 

in the Trotter direction), susceptibility and average energy of the system for various values of the 
transverse field at different temperatures. The antiferromagnetic order is seen to get immediately 
broken as soon as the thermal or quantum fluctuations are added. We discuss generally the phase 
diagram for the same LRTIAF model with perturbative Sherrington-Kirkpatrick (SK) type disorder. 
We find that while the antiferromagnetic order is immediately broken as one adds an infinitesimal 
transverse field or thermal fluctuation to the pure LRTIAF system, an infinitesimal SK spin glass 
disorder is enough to induce a stable glass order in the LRTIAF. This glass order eventually gets 

S destroyed as the thermal or quantum fluctuations are increased beyond their threshold values and 
the transition to para phase occurs. Analytical studies for the phase transitions are discussed in 
detail in each case. These transitions have been confirmed by applying classical and quantum Monte 
Carlo methods. We show here that the disordered LRTIAF has a surrogate incubation property of 
the SK spin glass phase. 
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Quantum phases in frustrated systems are being intensively investigated these days; in particular in the context of 
quantum spin glass and quantum axial next-nearest-neighbor Ising ( ANNNI) models [lH6( . Here we study in general 
the long-range Ising antiferromagnet put under transverse field (LRTIAF) with disorder in co-operative interactions 
superposed on it. We study here extensively, both analytically as well as numerically, the phase diagram for this model. 
As a special case, we also studied the pure long-range transverse Ising antiferromagnet model (i.e. no disorder). 

The finite temperature properties of sub-lattice decomposed version of the pure limit of this model was already 
considered earlier 0, [1] . The quantum phase transition and entanglement properties of the full long-range model at 
j_j zero temperature has also been studied @. 

Here we present some results obtained by applying analytical as well as Monte Carlo techniques [1(| to the general 
full long-range model at finite temperatures and transverse fields. We observe indications of a very unstable quantum 
antiferromagnetic (AF) phase (50% spin up, 50% spin down, without any sub-lattice structure) in the pure LRTIAF 
model, where the antiferromagnetically ordered phase gets destabilized by both infinitesimal thermal (classical) as 
well as quantum fluctuations (due to tunneling or transverse field) and the system becomes disordered or goes over 
to the para phase [llj. This kind of phase transition has also been studied by perturbative treatment (l2| . 

When a little spin glass-like disorder is incorporated with this pure LRTIAF model, the frustration is seen to 
destabilise the AF phase and stabilise a spin glass order. To check how this 'liquid'-like antiferromagnetic phase of 
the pure LRTIAF gets 'frozen' into spin-glass phase when a little disorder is added, we study in general the LRTIAF 
Hamiltonian with a coupling with the SK spin glass Hamiltonian and study this entire system's phase transition 
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behaviour induced by both thermal and tunneling field. Indeed, stable SK like spin glass phase is observed for both 
thermal or quantum fluctuations below finite threshold values [111 ]. 

This paper is organized in the following manner. In Section II, we introduce the quantum LRTIAF model with 
SK disorder. Then in the four subsections we discuss the analytical studies and simulations for some special cases 
of this general model. In Section II. A., we consider the quantum LRTIAF model without disorder, in Section II. B. 
the classical model (i.e., without the transverse field), in II. C. the quantum model at finite temperature and in II. D. 
the quantum model at zero temperature. The detail calculation of free-energy is given in Appendix A and an exact 
analysis at T=0 is given in Appendix B. In Section IV, we present some discussions on our results. 



II. LRTIAF WITH SK DISORDER 



The general model we study here is given by the following Hamiltonian 

1 N 

H = ~n E ( J ° + - hJ2*i - r£ of, (i) 

ij(j>i) i=l i 

where Jo is the parameter controlling the strength of the antiferromagnetic bias and J is an amplitude of the disorder 
Tij in each pair interaction, h and T denote respectively the longitudinal and transverse fields. The T controls the 
quantum-mechanical fluctuation. Here a x and a z denote the x and z component of the N Pauli spins 

< = (;_°i); <* = (!J) ; ' v 

As such the model has a fully frustrated (infinite-range or infinite dimensional) co-operative term. When we assume 
that the disorder Ty obeys a Gaussian with mean zero and variance unity, the new variable Jij = Jq + jTij follows 
the Gaussian distribution, P{Jij) — exp[— (Jy — Jo) 2 /2 J 2 ]/\/2~7r J. Therefore, we obtain the 'pure' antiferromagnetic 
Ising model with infinite range interactions when we consider the limit J — > keeping Jq < 0. Of course the model 
with Jo > and T = 0, is identical to the classical SK model and with Jo < and T = it is the LRIAF model. 

For an analytic (mean field) study of the model we define an effective magnetic field h e ff at each site, which is a 
resultant of the average cooperation enforcement in the z-direction and the applied transverse field in the x-direction, 
so that the above Hamiltonian can be written as, 



N 



H = h eff .J2^, (2) 



where 



i=i 



<Xj = of z + off, 



and 

heff = (h eff ) z z + (h eff ) x x 



h + J m z + J^/qy I + Fx, (3a) 



\heff\ = \J(h + J m z + JV^/) 2 +r 2 . (3b) 

This replacement of a z by its average value (a z ) = m z in (h e ff) z should be valid for this infinite range model (see 
Appendix A for its much more precise description under replica symmetric theory). The Gaussian distributed random 
field ^Jqy comes from the local field fluctuation (see e.g., [131 ]) given by the spin glass order parameter (see Appendix 
A for details). The average magnetisation is then given by 

Trae~ pH 



m = 



Tre-f 3H 



(tanh/31/^l).-^ (4) 

\ h eff\ 
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and hence we have 



£ 
£ 



By , Je " tanh/^/j 2 ,+T 2 (5) 
\/ J ^+ r2 

. r - tanh/^/j 2 , +T 2 (6) 

q = r By \ ~, Jeff \ tanh 2 pJj? ff + r 2 , (7) 

where J e // = ft + Jo mZ + J-J92/ an< ^ m = S»(°f) i s the magnetization and g = A -1 X^(°f ) 2 is the spin glass 
order parameter. We defined Dy = dy e~ y2 / 2 /V2ir. In Appendix A, we confirm that the above mean-field equations 
are identical to the results obtained by the replica symmetric theory at the ground state (/3 = oo). 

For the antiferromagnetic (Jo < 0) and/or the spin glass phase (with ft = 0), m z = is the only solution. We then 
have 

r 




tanh/3 A /(JVqy) 2 +r 2 , (8) 
/3\A J v%) 2 + r 2 . (9) 



II. A. Pure LRTIAF model 



The Hamiltonian of the infinite-range quantum Ising antiferromagnet (without any spin glass disorder) is 

N N 

i=\ i=i 

where Jo denotes the long-range antiferromagnetic (Jo < 0) exchange constant. We have denoted the co-operative 
term of H (including the external longitudinal field term) by and the transverse field part as H^ T K As such the 

model has a fully frustrated (infinite-range or infinite dimensional) co-operative term. 



H = + 

T N 

i,j(>i)=l 



II. A. 1. Analytical studies 



For h = 0, the Hamiltonian can be written as 

N N 

V 



H = - J i E 



N 

Without the disorder term, the h e /f in Eq. 3(a) gets modified to, 



i,j(>i)=l »=i 
/ N \ 2 AT 

# E-n - r E< • (id 



h e ff = Jom z z + Fx 



(12) 
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FIG. 1: Variation of the order parameter q (correlation in the Trotter direction) with transverse field r for T = 0.10,0.20 and 0.30 (h = 0) 
for two different system sizes (N = 100 and 200). q = for large I\ The inset shows the plot of q against the scaled variable T/T. 
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FIG. 2: Variation of the susceptibility \ with transverse field T for T = 0.10,0.20 and 0.30 (h < 0.1) for two different system sizes 
(N = 100 and 200). The corresponding susceptibility \cl f° r various temperatures for = 100 and 200 for the classical system are shown 
in the inset. \ converges to the classical values \cl f° r large T. 



With this modified field, the expressions of m z and m x become 



m> = (t^MJom'F + T*) V(Jo g + r2 d3a) 
m* = (tanh/gV^m^ + p) J . (13b) 



When Jo < 0, then m z — is the only solution of Eq. (13a). At zero temperature and at zero longitudinal and 
transverse fields, the would prefer the spins to orient in ±z directions only with zero net magnetization in the 

z-direction. This antiferromagnetically ordered state is completely frustrated and highly degenerate. Switching on 
the transverse field T would immediately induce all the spins to orient in the ^-direction (losing the degeneracy), 
corresponding to a maximum of the kinetic energy term and this discontinuous transition to the para phase occurs 
at T = (sec Appendix B for an exact result at T = 0). However, at any finite temperature the entropy term coming 
from the extreme degeneracy of the antiferromagnetically ordered state and the close-by excited states does not seem 
to induce a stability of this phase. 
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FIG. 3: Variation of average energy E with transverse field T for T = 0.10, 0.20 and 0.30 (h = 0) for two different values of N(= 100, 200). 
The corresponding average energy E c i for various temperatures for TV = 100 and 200 for the are shown in the inset. E converges to the 
classical values E c i for large T . 



II. A. 2. Monte Carlo simulation 



This Hamiltonian (fTO)) can be mapped to a (oo + l)-dimensional classical Hamiltonian [T(l[Ti| using the Suzuki- Trotter 
formula. The effective Hamiltonian can be written as (for J = — I) 



N M ^ N M 

j,j(>i)=l fe=l 2=1 fc=l 
N M 

+ ^ yiyi ^wi (14) 



i=i fc=i 

where 



B = (l/2)ln(coth(r/MT)). (15) 

Here M is the number of Trotter replicas and k denotes the A;-th row in the Trotter direction. B denotes the nearest- 
neighbor interaction strength along the Trotter direction. We have studied the system for N — 100. Because of the 
diverging growth of interaction B for very low values of Y and also for high values of M, and the consequent non- 
ergodicity (the system relaxes to different states for identical thermal and quantum parameters, due to frustrations, 
starting from different initial configurations), we have kept the value of M at a fixed value of 5. This choice of M value 
helped satisfying the ergodicity of the system up to very low values of the transverse field at the different temperatures 
considered T = 0.10,0.20 and 0.30. Starting from a random initial configuration (including all up or 50-50 up-down 
configurations) we follow the time variations of different quantities until they relax and study the various quantities 
after they relax. 

We studied results for three different temperatures T = 0.10,0.20 and 0.30 and all the results are for N = 100 and 
200 and M = 5. We estimated the following quantities after relaxation : 

(i) Correlation along Trotter direction (q) : We studied the variation of the order parameter 

N M 
i=l fe=l 

which is the first neighbor correlation along Trotter direction. Here, (...) indicate the average over initial spin 
configurations. This quantity q shows a smooth vanishing behavior. We consider this correlation q as the order 
parameter for the transition at r c . A larger transverse field is needed for the vanishing of the order parameter for 
larger temperature. The observed values (see Fig. [1} of Y c are ~ 1.6, 2.2 and 3.0 for T = 0.1, 0.2 and 0.3 respectively. 



G 



As shown in the inset, an unique data collapse occurs when q is plotted against T/T and one seems to get the complete 
disorder immediately as the scaling dos not involve any finite value T c . This is consistent with the observations in the 
next section. 

(ii) Susceptibility (\) : The longitudinal susceptibility x = (V-^-^O^Kj /c^fc)]/^' where h (— > 0) is the applied 
longitudinal held, has also been measured. We went up to h — 0.1 and estimated the x values. As we increase the 
value of the transverse held T from a suitably chosen low value, x initially starts with a value almost equal to unity 
and then gradually saturates at lower values (corresponding to the classical system where B = in Eq.(14)) as T 
is increased. Also at T = 0, the classical values are indicated in Fig. [2] This saturation value of x decreases with 
temperature. Again the field at which the susceptibility saturates are the same as for the vanishing of the order 
parameter for each temperature. 

(iii) Average energy (E) : We have measured the value of the co-operative energy for each Trotter index and then 
take its average E i.e. E = (H^ c ') of Eq. (10) with Jo = —1. It initially begins with —1.0 and after a sharp rise 
the average energy saturates, at large values of T, to values corresponding to the classical equilibrium energy (E c i 
for B = in Eq.(14)) at those temperatures. Again it takes larger values of T at higher temperatures to achieve the 
classical equilibrium energy. At T = 0, the corresponding classical values of E are plotted in Fig. [3J The variations 
of all these quantities indicate that the 'quantum order' disappears and the quantities reduce to their classical values 
(corresponding to B = 0) for large values of the transverse field T. 

The continuous transition-like behaviour seen from Fig. [I] can be justified from a mean field analysis (see Appendix 
B). At finite temperature it is the free energy that we have to minimise and the entropy term plays a crucial role. 
Minimization of free energy leads to an analytic variation of the total magnetization and no phase transition at any 
finite temperature. 

II. B. LRIAF with disorder : Classical model at finite temperature 
II.B.l. Analytical studies 

For the classical case, i.e., T = 0, Eqns. §5$ and |7]) reduce to (for h = 0) 



m = I Dy tanh f3 (J \fqy + J m) (17) 

J — OO 

q = / L>2/tanh 2 (/3JV^ + J m) , (18) 



where m z has been replaced by m. For Jo < 0, again we find that m = is only physical solution for all temperature 
regimes. This means that there are three possible phases : namely, the antiferromagnetic phase, the paramagnetic 
phase and the spin glass phase. In all these three phases, the magnetization m is zero. To determine the critical point 
Tsg at which the spin glass transition takes place, we expand the equation with respect to q for q ~ and m — 0. In 
the limit of J ^ 1, we have 

,= """>'- 1 . (19) 
2(/3J) 4 

We therefore have Tsg = J and the critical point is independent of the antiferromagnetic bias Jo. This result means 
that the antiferromagnetic order can appear if and only if we set Jo < and T/J = 0, Jo/ J = — oo. On the other 
hand, for — oo < J < at low temperature regime T < Tsg, the spin glass phase appears. We plot the phase diagram 
in Fig. [4] We also conclude that the system described by the Hamiltonian (fit)]) with T = is immediately 'frozen' 
when we add any infinitesimal disorder J > 0. 

From the view point of the degeneracy of the spin configurations, we easily estimate the number of solution for the 
antiferromagnetic phase as An/[(7V/2)!(7V/2)!] ~ 20.693^ ^ n ^ e antiferromagnetic ground state only N/2 spins out 
of N will have free choice (to be up or down) and the rest have to follow), which is larger than the number of the 
SK model e 0199Ar [l3j |. However, for the infinite range antiferromagnetic model, the energy barrier between arbitrary 
configurations which gives the same lowest energy states is of order 1 and there is no ergodicity breaking. 
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FIG. 4: The phase diagram of classical SK model |l3j extended for antiferromagnetic bias. For Jo < 0, there exist spin glass phase 
below T/J = 1 and the critical temperature is independent of the strength of the antiferromagnetic bias Jo. For pure LRIAF at finite 
temperature (T > 0), the anti- ferromagnetic order disappears and the system changes to the paramagnetic phase. When we add an 
infinitesimal disorder J > 0, the antiferromagnetic order is broken down and the system suddenly gets 'frozen' into a spin glass (SG) 
phase. 



II.B.2 Monte Carlo studies 



In order to study the effect of introducing disorder in the classical LRIAF model, we performed a Monte Carlo study 
with a system of N — 100 spins. But the distribution function (P) of disorder introduced here is different. Instead of 
a Gaussian distribution, we applied a binary distribution (± J) with a probability p : 

P(Jij) = p5{J ij - J) + (1 - p)S{J ij + J) . (20) 

In our study we kept J — 1 and Jo = J, where J is the average interaction strength. Each of the Ising spins interact 
with every other ferromagnetically with probability p and antiferromagnetically with probability 1 — p. Therefore 
J = 2p — 1 and J — 2y / p(l — p). At the two limiting values of p — and p = 1, the system becomes purely 
antiferromagnetic (LRIAF) and purely ferromagnetic respectively. Thus at these two limiting values of p, the system 
has no fluctuation at all (J = 0). Whereas for p = 0.5, the fluctuation is maximum (J = I). So as predicted above, 
we need maximum thermal fluctuation to destroy the glassy phase. 

To identify the glass phase, we considered a replica of the original system to be studied and evolved the two systems 
simultaneously by Monte Carlo technique. We also repeated the study for N = 50 and 200 also. For all sizes we 
obtained almost same phase diagram. We measured the average absolute value of the spin-spin correlations of the two 
systems (the original and the replica one) at different times for a certain temperature. Let us denote this parameter 
by q and a and f3 denotes the original system and the replica. 



1 N 

— Vs a / 
N ^ 1 1 

i=l 



(21) 



We measure the average steady state values of this parameter for various temperatures. The average is over different 
configurations. When p is very high i.e., when the number of ferromagnetic bonds are dominant, the ferromagnetic 
to paramagnetic transition can be identified easily by measuring the change in magnetisation and divergence of 
susceptibility at the critical temperature. But for p < 0.5 and slightly above 0.5, the magnetisation remains low for all 
temperatures and hence it is difficult to identify the existence of any glass phase/paramagnetic phase seperately. The 
glass and paramagnetic phase can be distinguished by studying the steady state values and fluctuation of q. We have 
studied for three different values of p — 0.20,0.30 and 0.40. For a particular value of p, upto a certain temperature 
the value of q remains high indicating the spin glass phase. The fluctuations in q remain small. This implies that the 
flipping of spins is very low (freezing of spins). As the temperature is raised, the steady state value of q decreases 
gradually (inset of Fig. [SJ but the fluctuation increases. At a certain temperature (glass-para phase transition point) 
this fluctuation reaches maximum and above this temperature the value of q gradually goes to zero. With further 
increase of temperature, the fluctuation decreases indicating a second order glass-para phase transition (Fig. [5]). 

It can be observed from (Fig. [5]), that as we increase the value of p from to 0.5, the transition temperature 
increases. It is expected from our previous explanation, as p approaches 0.5 ^fluctuations J increases and so is the 
transition temperature. We have given a plot of the T c /J as a function of J/ J (Fig. 
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FIG. 5: Variation of the fluctuation of spin-spin correlation q with temperature T for p = 0.40, 0.30 and 0.20 for N = 100 for classical 
disordered system. The inset shows the plot of q against temperature T. 
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FIG. 6: The phase diagram of classical LRIAF with disorder obtained from Monte Carlo simulation. This phase diagram is similar to 
that obtained from static and replica symmetric approximations. 

II. C. LRTIAF with disorder : Quantum model at finite temperature 
II.C.l. Analytical studies 

The approximate saddle point equations have already been presented in equations [5j [6] and [7] For detail calculations 
see Appendix A. The variations of m x , q and q are shown in Fig. [7] The phase boundary between the spin glass and 
paramagnetic phases is given by setting m z = and q ~ and we get 



Jtanh 



T 



(22) 



Obviously, the boundary at T = gives Tsg = J- On the other hand, when we consider the case of T ~ 0, we have 
Tsg = J (consistent with the classical result). These facts imply that there is neither an antiferromagnetic nor a spin 
glass phase when we consider the pure case J = because the critical point leads to Tsg = ^sg = 0. Therefore, we 
conclude that the antiferromagnetic phase can exist if and only if T = T = (Fig. |S| . 
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FIG. 7: The result of numerical calculations for the saddle point equations for m x , q and q as a function of T for T^O. 
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FIG. 8: Phase diagram for the quantum system. The antiferromagnetic order exists if and only if we set T = V = and J = 0. As 
the J decreases, the spin glass phase gradually shrinks to zero and eventually ends up at an antiferromagnetic phase at its vertex (for 
r = = T = J) as discussed in Section II. C.I. 

II. C. 2. Quantum Monte Carlo studies 

Now to study the quantum system we again apply the finite temperature quantum Monte Carlo method as applied 
to study the pure LRIAF model in Section II A. Like the classical system, here also we incorporated a disorder with 
binary distribution (±J) with a probability p. As mentioned earlier, to study by quantum Monte Carlo, we map 
the Hamiltonian (|T|) to a (oo + 1) dimensional classical Hamiltonian using the Suzuki- Trotter formula. The effective 
Hamiltonian can be written as (for h — 0) 

N M N M 

^ = _ JVM ^ /J JijVi,kGj,k + B ///, Q"»,fcQ-j,fc+i , (23) 

i,j(>i) = l fc=l i=l fe=l 

where 

B = (l/2)ln(coth(r/MT)). (24) 

M, k and B represents the number of Trotter replicas, fc-th row in Trotter direction and nearest-neighbour interaction 
strength along the Trotter direction respectively. The distribution P(Jy) is given by Eq. (20). To identify the 
glass phase, a replica of the original system to be studied has been taken and the two systems has been evolved 
simultaneously by Monte Carlo technique. 
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FIG. 9: |(a)| ^ariation of the fluctuation of spin-spin correlation q w ith V at T = 0.1 for p = 0.40, 0.30 and 0.20 for N = 100 
(M = 5). The inset shows the plot of q against transverse field I\ (b) Variation of the fluctuation of spin-spin correlation q 
with r at T = 0.2 for p = 0.40, 0.30 and 0.20 for TV = 100. The inset shows the plot of q against transverse field T. 
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FIG. 10: Variation of r c with J for the quantum LRIAF with disorder obtained from Monte Carlo simulation. This plot is qualitatively 
similar to that obtained from static and replica symmetric approximations. 



The quantity measured (q) is the average absolute value of the spin-spin correlations of the two replica systems at 
different times for a certain temperature and T, where a and /3 denotes the original one and the replica. 



1 



NM ^ 

i=i 



JVM 



(25) 



Initially we kept the temperature fixed at a certain value (typically T = 0.1) and measured the average steady state 
values of this parameter for various values of T. The average is over different configurations. Here also the glass and 
paramagnetic phase can be distinguished by studying the steady state values and fluctuation of q. For a particular 
value of p, upto a certain T the value of q remains high indicating the spin glass phase. The fluctuations in q remain 
small. This implies that the flipping of spins is very low (freezing of spins). As the value of T is raised, the steady 
state value of q decreases gradually (inset of Fig. 9(a) ), but the fluctuation increases (Fig. 9(a) I. At a certain value of 
r (say T c ) (glass-para phase transition point) this fluctuation reaches maximum and above this T c the value of q goes 
to zero. With further increase of T, the fluctuation decreases indicating a second order glass-para phase transition. 
The observed values of r c are ~ 0.2, 0.35 and 0.58 for p = 0.2, 0.3 and 0.4 respectively. A similar study for T = 0.2 is 
presented in Fig. |9(b)[ The values of T c are ~ 0.10,0.20 and 0.40 for p = 0.2,0.3 and 0.38 respectively. 

It has been noticed, that as we increase the value of p from to 0.5, the transition field increases. It is expected 
from our previous explanation, as p approaches 0.5, fluctuations J increases and so is the transition field T c . We have 
given a plot of the r c as a function of J (Fig. [TU]) . Though these results do not satisfy the predictions of static and 
replica symmetric approximations quantitatively, but qualitatively they reflect all the features. 
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IID. LRTIAF with disorder : Quantum model at zero temperature 



We study this case only analytically. As well-known, in the mean-field description of the pure quantum transverse 
Ising systems, the total magnetization is conserved as (m z ) 2 + (m x ) 2 — 1. However, if some disorders are taken into 
account, it is non-trivial problem to answer the question; if the magnetization conservation still holds or not. As we 
discussed before, for antiferronmagnets, m z is always zero and the magnetization conservation reads (m 1 ) 2 = 1. In 
following, we derive the condition on which the magnetization conservation holds. 

For finite transverse field but zero temperature, i.e., /3 = oo, equations (8) and (9) reduce to 

By -, F = (26) 

Du{ , JV ~ qV } (27! 

which are obtained within the replica symmetric theory in this zero temperature limit (see Appendix A) . The variations 
of m x and q are shown in Fig. [TT] It can also be seen that, for no disorder, i.e. J = 0, we obtain m x — 1 and q = 0. 




FIG. 11: The result of numerical calculations for the saddle point equations for m x and q as a function of T for T = 0. 

For finite disorder J ^ 0, we expand equation ([27| with respect to J^/qy <C 1. Then, we have 1 ~ (J/r) 2 — 3q( J/r) 4 , 
namely, 




where we used £° Dyy 2 = 1, £° Dyy A = 3. Eqn. ^ implies q ^ for T < J. Hence, 




From this result, we find that magnetization conservation is broken for T < J. This fact means that each spin starts 
'frozen' locally and the degree of freedom of spins is remarkably reduced. Therefore, m x might be a good indicator to 
detect the spin glass phase for the system in which the both order (antiferromagnet) and disorder (spin glass) phases 
possess the same spatial symmetry with m z = 0. 
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FIG. 12: Variation of m x and q with J/T 



III. DISCUSSION 



We considered here first a long-range Ising antiferromagnet with disorder incorporated in it and put it in a transverse 
field. Although long-range interaction is unrealistic and also difficult for numerical studies, it is convenient for 
analytical approaches such as mean-field calculations which are exact in certain limits. So here we have studied this 
kind of long-range interaction. We have obtained the finte temperature free-energy expression (Appendix A) for this 
model and studied analytically the magnetisation, spin glass order and the correlation (Trotter replica). For the 
pure case (i.e. no disorder) the antifcrromagnetic order is seen to get immediately broken as soon as the thermal or 
quantum fluctuations are added (see Sec. II. A.). However, when we add the disorder as in SK Hamiltonian, on that 
of LRTIAF as perturbation, we find that an infinitesimal spin glass disorder is enough to induce a stable glass order 
in this LRTIAF antiferromagnet (Sec. II.B.). This glass order eventually gets destroyed as the thermal or quantum 
fluctuations increased beyond their threshold values and the transition to para phase occurs (Sec. II.B. and II. C). As 
shown in the phase diagram in Fig. SI the antiferromagnetic phase of the LRTIAF (occurring only at J = = T = T) , 
can get 'frozen' into spin-glass phase if a little SK-type disorder is added (J ^ 0); the only missing element in the 
LRTIAF (which is fully frustrated, but lacks disorder) to induce stable order (freezing of random spin orientations) 
in it. These results has been confirmed by Monte Carlo simulations. We have not addressed the question of possible 
replica symmetry breaking in our study. The degeneracy factor e°- 693N of the ground state of the LRIAF is much 
larger than that e 0199N for the SK model. Hence, (because of the presence of full frustration) the LRIAF possesses 
a surrogate incubation property of stable spin glass phase in it when induced by addition of a small disorder. 
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In this appendix, we show the derivation of the free energy per spin for the system to be described by the Hamil- 
tonian. 
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Appendix A: Derivation of free energy 
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Carrying out the Suzuki- Trotter decomposition, we have the replicated partition function. 



z m = tr w exp 



B = - In coth 
2 



Tf E E E '«°f (*K(*) + B EEE °f (*K(* + 1) 



M 
M 



i k a 



(A2) 



(A3) 



where a and fc denote the replica and Trotter indices. M is the number of the Trotter slices and f3 is the inverse 
temperature. The disorder obeys 



1 



cxp 



(Jjj-jo) 2 
2 J 



(A4) 



In other words, the follows 



Jij = jo + Jx, P(x) = —=e x 2 . (A5) 

V 2vr 

We should notice that jo > 0, J = is pure ferromagnetic transverse Ising model, whereas jo < 0, J = corresponds 
to pure antiferromagnetic transverse Ising model. Then, by using Dxe ax = e a / 2 , Dx = dxcT x / 2 /\/~2n, we have 
the average of the replicated partition function as 



< ZIj > = tr {CT} exp 



M 



E E E ^?( fc K( fc ) + B EEE + 1) 



fc a ij 



fc a i 



x exp 



/3 2 J 2 
2M 2 



E EE°?(*KW 



tr {fT} exp 
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E E E ^?( fc K( fc ) + B EEE *?(*K(* + 1) 
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x exp 



/3 2 J 2 
2M 2 



(A6) 



where the bracket was defined as <C • • • >•= JYlij dJijP(Jij)(- ■ ■ ). To take a proper thermodynamic limit, we use 
the scaling 



Jo . J 



(A7) 



For this rescaling of the parameters, the averaged replicated partition function <C ^> reads 



«^» = tr W / nn 



dm a (k) 



k a 



y/2nM/pj N 
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x exp 
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2M 
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We next assume the replica symmetry and static approximations such as 

m a (k) = ( ( r tt (fc)>-^E^( fc ) = m 



N ■ 

q aP (k,k') = (a a (fc)^(fc')) = l^af(fc)af(fc')=g 

i 

q aa (k,k) = (a a (k)a a (k)) = 1 £ of (fc)af (fc') = 9 



Then, we should notice the relation: 

_\ 2 



«EEE°f(*)°f(*') - (#) C(EE^)) «Ee(e<w 
ffEEE^(*K(*') - te)^EE(E<( fc )) 2 

fcfc' a i \ / » a \ fe / 



To take into account the above relations, we obtain in the limit of AT — » oo as 

dm a (k) Z" 00 tttt dq a p(k,k ) f°° yryj dq aa (k,k ) 



«^m»= nn^=== nn nn 



x exp 



/3Jo 2 (/3J) 2 2 (/3J) 2 , 2 
niV — m H — q — q 



2 4 

exp[nNf]. 



/■OO pOO i 

Dyln / Du 2 cosh/3^/ ( J m + J^y + J^q - qu) 2 + V 2 
-oo J — oo 

Therefore, the following / is regarded as free energy per spin by the definition of replica theory 

f = i 

J 2 4 y 4 1 

/OO f>00 , ■ 

Dyln / Du 2 cosh/3^/ ( J m + J^y + Jy/q - qu) 2 + T 2 
-oo J —oo 



1. Saddle point equations 

For simplicity, we define 



b = J m + J^fqy + JyJ q- qu 
9 = Vb 2 + r 2 



Then, we have the following simplified free energy 

Dy\n [ 

J —OO J - 



f = -^ m 2 + i^Dl( q 2 -q 2 )+ I Dyln / i? U 2cosh/3e. 
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The saddle point equations are derived as follows [H, [H| . 
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2. At the ground state 

We hrst should notice that q is always larger than q. In fact, we can easily show that 

7.!°^ (|r) cosh/36^ 
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Then, we consider the limit of /? — > 00. If q — q = e > is of order 1 object, the free energy / diverges in the limit 
of /3 — > 00 as (/3J) 2 (q 2 — g 2 )/4. Therefore, we conclude that q = q should be satisfied in the limit of /3 — > 00 and we 
obtain the saddle pint equation at the ground state as 



m T = 



Dy [e 



Dy 



Dy 
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(J m + Jy/qy) 
(Jom + Jy/qy)* + T 2 

(J m + Jy/qy) 



Dy 



Dy 



l{J m + J^qy) 2 + T 2 
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(J m + J^qy) 2 + r 2 
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Appendix B: Exact analysis at T = 

First, let us consider the case of pure LRIAF model and rewrite our Hamiltonian H in Eq.(l) for J = as 



H 



/ N \ 2 N N N 

2^ E-n -^E«) 2 -»E<- r E^ 

Vi=l / i=l i=l i=l 



(Bl) 



If we now denote the total spin by a to t i-e. Stot — js X^:=i &i (where N\a\ = 0, 1, 2, N), then the Hamiltonian H 
can be expressed as 



H 1 2 



tot 1 °tot 



1 



(B2) 
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Let us assume the average total spin (a) to be oriented at an angle 9 with the z-direction : (of ot ) = mcos9 and 
(af ot ) = m sin 9. Hence the average total energy E to t = (H) can be written as 

— ^ = \m 2 cos 2 9 — hmcos9 — Tmsin9 — (B3) 
N 2 N K J 

At the zero temperature and at T — 0, for h = 0, the energy Etot is minimised when = and m = (complete 
antiferromagnetic order in z-direction). As soon as T ^ (h = 0) the minimisation of Etot requires 9 = n/2 and 
m = 1 (the maximum possible value); driving the system to paramagnetic phase. This discontinuous transition at 
T = was also seen in [9(. As observed in our Monte Carlo study in the previous section, T C (T) — > as T — > 0. 
This is consistent with this exact result r c = at T — 0. For T = (and h = 0), therefore, the transition from 
antiferromagnetic (6 = = m) to para (0 = 7r/2, m = 1) phase, driven by the transverse field T, occurs at T = itself. 

One can also calculate the susceptibility \ a t T = = T. Here E to t/N = ^m 2 cos 2 9 — hmcos9 — -4 and the 
minimisation of this energy gives mcos9 = h giving the (longitudinal) susceptibility \ — mcos9/h — 1. This is 
consistent with the observed behaviour of x shown in Fig. [2] where the extrapolated value of x at T = increases 
with decreasing T and approaches x = 1 as T — > 0. 

At finite temperatures T ^ 0, for ft, = 0, we have to consider also the entropy term and minimise the free energy 
J- = E to t — TS rather than E to t where 5* denotes the entropy of the state. This entropy term will also take part in 
fixing the value of 9 and m at which the free energy T is minimised. As soon as the temperature T becomes non-zero, 
the extensive entropy of the system for antiferromagnetically ordered state with m ~ (around and close-by excited 
states with 9 = 0) helps stabilisation near 9 = and m = rather than near the para phase with 9 = -k/2 and 
m = 1, where the entropy drops to zero. While the transverse field tends to align the spins along x direction (inducing 
9 = Tt/2 and m = 1), the entropy factor prohibits that and the system adjusts 9 and m values accordingly and they 
do not take the disordered or para state values (9 = n/2 and m = 1) for any non-zero value of T (like at T = 0). 
For very large values of T, of course, the free energy T is practically dominated by the transverse field term in H 
and again 9 = 7r/2 and m = 1, beyond T = T C (T) > for T > 0. However, this continuous transition- like behaviour 
may be argued [17| to correspond to a crossover type property of the model at finite temperatures (suggesting that 
the observed finite values of r c (T) are only effective numerical values). In fact, for h = one adds the entropy term 
— T\nD, where D is the degeneracy for getting total spin r = N\3tot\ @, 

Nl Nl 
D ~ (N/2 + t)1(N/2-t)1 ~ (N/2 + t + 1)\(N/2-t- 1)1 ^ 

to E^t in Eq.(A3) to get T and one can then get, after minimising the T with respect to m and 9, m = tanh(r/2T), 
which indicates an analytic variation of m and no phase transition at any finite temperature for J = (antiferromag- 
netic phase occurs only at T = T = and Jq/J = oo as shown in Fig. @|. 
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